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Abstract 

In this paper we have examined the internal dynamics of the globular cluster NGC 362 
using a combination of V band CCD images and echelle spectra of member red giants. A 
V band surface brightness profile (SBP) was constructed from the CCD images, and, after 
it was determined that the cluster is not post core-collapse, fit with single- and multi- 
mass King-Michie (KM) models. We found that for small values of the mass function 
slope, x, anisotropic models were favored while for steeper mass functions isotropic orbits 
provided superior fits. The total cluster luminosity is 1.70 ±0.1 x 10 5 Ly© [assumes (m- 
M)o = 14.77]. A total of 285 stellar spectra were obtained of 215 stars for radial velocity 
determinations. Three stars were obvious non-members and four showed strong evidence 
for radial velocity variations; these latter stars are probably members of binary systems 
with periods less than a few years. The true cluster binary fraction was determined from 
simulations to be 0.15 for circular orbits or 0.27 for orbits with a distribution function 
/(e) = e (e is eccentricity). This relatively high binary detection frequency may indicate 
that NGC 362 is overabundant in binaries compared to other clusters. The 208 remaining 
stars showed no sign of rotation and had kinematics which were incompatible with KM 
models having isotropic orbits and luminosity profiles consistent with the SBP. Therefore, 
the best agreement with both the kinematic data and the SBP were for shallow mass 
functions x = 0.0 — 0.5 and intermediate amounts of anisotropy in the velocity dispersion 
tensor. In this best-fit range, the derived cluster mass is M = 2.5 — 3.5 x 10 5 M© for a global 
mass-to-light ratio of M/Ly = 1.5 - 2.0 M©/Ly©. This low value for x is in disagreement 
with the correlation between x and the height above the Galactic disk seen for a sample 
of other clusters. The results are also different from the sharp turn-up in the low mass 
end of mass functions derived from some deep luminosity functions of three other globular 
clusters. 
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1. Introduction 

Kinematic studies are a potentially powerful means of learning about the present 
dynamical state and stellar content of globular clusters. They can lead to important 
conclusions regarding the history of the globular cluster system and the low-end of the 
stellar mass function which impact on star formation and the formation and evolution of 
the Galactic halo. 

The exact nature of the low mass end of the halo mass function is highly uncertain. 
Obtaining this information through observations of individual field halo stars is very diffi- 
cult as one is faced with contamination from faint galaxies, uncertain distances and small 
number statistics. One can overcome these problems by studying the luminosity function 
of halo stars which are members of globular clusters. Unfortunately these studies are gen- 
erally confined to small regions at large projected radii from the cluster centers and, with 
a very few noteworthy exceptions, sample only a small range in mass (i.e. 0.5 - 0.8 M©). 
Other problems include crowding, background contamination and corrections to the local 
mass function for mass segregation. 

Kinematic measurements along with dynamical modeling can be used to overcome 
some of the problems in the luminosity function approach, but also possess some intrinsic 
shortcomings. The Achilles heel of this method is the problem of uniqueness which occurs 
whenever one is dealing with a system where mass does not follow light (galaxy, and 
galaxy clusters share this problem). In such cases, it is possible to obtain a large range of 
models possessing different dynamical parameters (i.e. mass, orbital anisotropy) which are 
similarly consistent with the observational data. The approach that one makes, therefore, 
is to restrict the parameter space to a "reasonable" subset of the possible range. 

Until recently, obtaining large samples of high precision radial velocities has been a 
very time-consuming task. Thus, there are very few kinematically well-studied clusters 
(i.e. possessing more than 100 or so radial velocities). Among the better-sampled clusters 
are M3 (Gunn & Griffin 1979), Omega Cen and 47 Tuc (Meylan & Mayor 1986), M13 
(Lupton, Gunn, & Griffin 1987), M15 (Peterson et al. 1989), NGC 6397 (Meylan & Mayor 
1991), NGC 3201 (Cote et al. 1993), with several large data sets in existence yet to be 
published. 

In this paper we present a dynamical analysis of NGC 362. This cluster is an ideal 
candidate for a radial velocity study because of its high systemic radial velocity (223.5 km 
s _1 ) which enables one to unambiguously remove contaminating field stars. Furthermore, 
it is a particularly difficult cluster for which to obtain a luminosity function because it is 
projected on to the halo of the Small Magellanic cloud. In §2.1 and 2.2 we discuss the 
CCD imaging and echelle spectroscopy observations and reductions. In §3.1 and 3.2 we 
describe the "reasonable" dynamical parameter subspace we have chosen for this analysis, 
that is the popular King-Michie models which invoke mass segregation, tidal truncation 
and velocity anisotropy with a minimum of free parameters. In §3.3 we summarize the 
results of the modeling, including a comparison with some previous luminosity and mass 
function studies. §3.4 is a justification of the use of the King-Michie models based on the 
NGC 362 relaxation timescales. 



2. Observations and Reductions 
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2.1 Surface Photometry 

Four V CCD frames of NGC 362 were obtained at the Las Campanas Observatory 
(LCO) 1.0 m telescope on 1991 January 22. The TEK2 1024 2 chip was used (readout noise 
= 7 e~, gain = 2 e~/ADU, and angular scale = 0.61" px _1 ). The exposure times were 30 
seconds for each image. The frames were positioned with the cluster in one corner of each 
image such that full azimuthal coverage was obtained. 

The coordinates of the cluster center were adopted from Shawl & White (1986). As- 
trometric zero-points for each frame were determined through the cross-identification of 
15 stars from the Guide Star Catalog (Lasker et al. 1990, Russell et al. 1990, and Jenkner 
et al. 1990). The RMS uncertainty in the position of a selected Guide Star was 0.5 arcsec. 

Surface photometry was performed in a manner similar to Djorgovski (1988). The 
frames were broken up into a series of concentric circular annuli centered on the cluster. 
The annuli were further divided into eight azimuthal sectors. The average pixel brightness 
was determined for each sector in a given annulus and the median of the eight separate 
measurements was taken as the representative brightness at the area-weighted average 
projected radius of the annulus (i.e., the mean radius of all the pixels within the annulus 
which is approximately equal to the geometric mean). The standard error of the median 
of the eight sectors was adopted as the photometric uncertainty. 

A background level (a combination of SMC stars, sky light and Galactic foreground 
stars) was estimated from regions at large projected distances from the cluster. We found 
that the surface brightness profile (SBP) tended to level out beyond 21 pc (we have adopted 
a distance to the cluster of 9 kpc, see below). By "leveling out" we don't necessarily mean 
that the cluster light does not extend beyond this point but simply that fluctuations in the 
background dominate to such an extent that it is no longer possible to observe the profile 
declining in intensity. Therefore, it was this region with a projected radius of 21.0 < R 
(pc) < 33.0 that was used for the background determinations. The background value was 
195.7 ± 1.6 L y0 pc" 2 . 

The background-subtracted surface photometry data is presented in Table 1 [assuming 
a cluster distance of 9 kpc (see below), My = 4.83 Mihalas and Binney 1981, p. 60]. 
Column 1 is the area-weighted angular radius, column 2 is the projected area-weighted 
radius, column 3 is the V luminosity surface density, and columns 4 and 5 are the inner 
and outer radii of the annuli, respectively. Calibration was accomplished by matching up 
10 on-frame standards from Alcaino et al. (1988). The RMS scatter was 0.055 mag. 

2.2 Radial Velocities 

Spectra of red giants and horizontal branch stars in the region surrounding NGC 362 
were obtained during two runs: 1989 December 7 - 14, and 1991 December 15 and 17, 
using the photon-counting echelle spectrograph on the 2.5m Dupont reflector, designed and 
built by Steve Shectman. The observation and reduction procedures have been discussed 
extensively in Welch et al. (1991) and remain largely unchanged for these data. Briefly, the 
observing procedure consisted of exposures with integration times of 100 - 500s and Th-Ar 
arcs approximately every 45 minutes. A representative LCO spectrum is shown in Fig. 2 
of Cote et al. (1991). The reduction utilizes the IRAF ECHELLE and RV packages (Tody 
1986) to obtain both velocities and velocity uncertainties according to the prescription of 
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Tonry & Davis (1979). The velocity zero-point is tied to the IAU velocity standard 33 Sex 
as described in Fischer et al. 1992 and is accurate to better than 2 km s _1 . 

A total of 285 stellar spectra were obtained. There were 215 distinct stars observed; 
fifty were observed twice, one observed three times, one observed ten times, and one was 
observed eleven times. The heliocentric radial velocities are listed in Table 2, where the 
columns from left to right are: the star designation, the equinox J2000.0 right ascension 
and declination, the derived projected radius in arcsec and position angle in degrees, the 
heliocentric Julian Date, the observed radial velocity and uncertainty, the mean radial 
velocity and uncertainty for stars with multiple observations, the x 2 an d number of degrees 
of freedom for stars with multiple observations, and comments. The designations are taken 
from Harris (1982) when possible. The adopted absolute distance modulus (see §3.1) of 
(m-M)o = 14.77 mag corresponds to 0.0436 pc/arcsec. 

As pointed out above, our complete data set includes 53 stars having radial velocities 
corresponding to two different epochs (separated by roughly two years). Such a timespan 
should be sufficient for identifying cluster binaries, and, indeed, four of these 53 stars (or 
7.5%) show velocity differences greater than 6 km s -1 . Motions in the atmospheres of stars 
near the tip of the RGB rarely exceed about 5 km s _1 (Mayor et al. 1984). We therefore 
suspect these discrepancies arise from the effects of a companion since an inspection of the 
catalog of (known) NGC 362 variables (Sawyer-Hogg 1973) revealed only one such star in 
our sample (H1204=V11). If we insist on restricting ourselves only to those stars showing 
variations greater than 8 km s -1 , our number of detected binaries drops to three (or 5.7%). 
In either case, the NGC 362 binary detection frequency appears significantly higher than 
that seen in other globular clusters (Hut et al. 1992). 

In order to convert these detection frequencies into true binary fractions, we have 
carried out a series of Monte Carlo simulations similar to those of Pryor et al. (1988). 
We have generated 10000 simulated radial velocities pairs (separated by two years) for a 
grid of hypothetical clusters having binary fractions, Xb = 0.00, 0.05, 0.10, 0.15, 0.20, 0.25, 
0.30, 0.35, 0.40, 0.45 and 0.50. Unlike the models of Pryor et al. (1988), which include a 
second free parameter, r m i n , (the minimum orbital separation which avoids mass transfer) 
we have randomly assigned an r min to each binary by choosing values in the range 0.05 
- 0.14 AU. The vast majority of our program objects uniformly lie in the range 14.0 < 
V < 15.5 which, according to the [Fe/H] = -1.03, [O/Fe] = 0.5, T = 13 Gyr isochrone 
of Bergbusch and VandenBerg (1992) (see §3.1), corresponds to the above primary radii. 
Simulations have been carried out for the two cases of (1) purely circular orbits and (2) a 
distribution of orbital eccentricities given by /(e) = e (Heggie 1975). The models use the 
distributions of orbital periods and mass ratios described by Pryor et al. (1988). In all 
cases, we assume random orbital inclinations, orientations, phases and observational error. 

In Table 3 we show the fraction of stars expected to show velocity variations greater 
than 6 or 8 km s _1 based on these simulations. The first column gives the assumed binary 
fraction while the second and third columns record the expected fractions for both circular 
and eccentric orbits, respectively. By comparing our observed detection frequency to the 
model results, we see that the true NGC 362 binary fraction is roughly 0.15 for circular 
orbits or 0.27 for a Heggie (1975) distribution of orbital eccentricities. Of course, these 
numbers may be slightly enhanced by our procedure of preferentially selecting for repeat 
measurement those stars which showed higher than normal residual velocities. Neverthe- 
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less, the derived binary frequency for NGC 362 seems rather higher than that seen by Hut 
et al. (1992). Clearly, further monitoring of these four stars (H1348, H1419, H2205 and 
H2222) is in order, especially the luminous RGB star H2205 which shows extreme velocity 
variability (but no known photometric variability whatsoever (Alcaino 1976 and Harris 
1982). 

For the following dynamical analysis, we have therefore rejected seven stars from 
the sample in Table 2. In addition to the four binaries mentioned above, the RR Lyrae 
H1204 and the two (photometric and radial velocity) non-members H1423 and H2113 
have also been discarded. For the remaining stars which have repeat measurements, the 
total x 2 is then 70.47 for 67 degrees of freedom which corresponds to P(x 2 ) ~ 0.35. The 
analytical error estimates returned by RVXCOR therefore seem to be reasonable (for this 
particular data set anyway) and we do not need to add a so-called "jitter" term which other 
investigators have found necessary to account for possible motions in the atmospheres of 
these evolved stars. It is possible that our error estimates are more generous and therefore 
take into account this atmospheric jitter or simply that other studies, which have generally 
relied on radial velocity scanners, have tended to underestimate the velocity uncertainties. 

The 208 remaining radial velocities are plotted in Fig. 1 versus projected radius (top 
panel) and position angle (bottom panel). There is no obvious indication of rotation in 
the velocity versus position angle graph which should manifest itself as a sinusoidal trend 
if aligned close to equator-on. An alternative way to test for rotation is to measure the 
difference in median velocities for stars on either side of an imaginary axis as it is stepped 
around the cluster center in l°intervals. Fig. 2 is a plot of this and the solid line is the 
best-fitting sine curve having an amplitude of 0.38 km s _1 . While there does seem to be 
a correlation in the data, we found, based on 1000 simulations of the data assuming a 
non-rotating model, that this amplitude is not significant. That is, one would expect this 
amplitude or higher over 90% of the time if the cluster was not rotating. Therefore, we 
conclude that if rotation is present in NGC 362, it is too small an effect compared to the 
velocity dispersion for us to detect (see §2.2). Furthermore, if rotation is present, it is 
not dynamically important, meaning that the assumption of non-rotation and sphericity, 
consistent with the small cluster ellipticity of e = 0.99 found by White & Shawl (1987), 
will not significantly bias our mass estimates. If the cluster is being viewed pole-on or at 
a significant inclination, rotation would not be observable even if present. However, the 
small amounts of flattening seen in Milky Way globular clusters indicates that rotation is 
probably dynamically insignificant in most of them. 

Finally, Suntzeff et al. (1986) reported lower precision velocities for four stars in NGC 
362: H1216 (242.3 km s" 1 ), H1441 (210.8), H2108 (225.8), and H2431 (220.1). Our ve- 
locities for these stars are different by —15.0, +18.5, —1.5, and +4.6 km s , respectively, 
which suggests good agreement in the velocity zeropoint considering that their estimated 
uncertainty for an individual velocity was 15 km s _1 . 

3. King-Michie Models 

We used maximum likelihood fits of single- and multi-mass King-Michie (KM) models 
(King 1966, Michie 1963, Da Costa & Freeman 1976 Gunn & Griffin 1979) to the SBP in 
order to constrain the orbital distribution function and obtain the cluster luminosity. Once 
this is accomplished, maximum likelihood scaling is determined from the radial velocities 
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yielding cluster mass estimates. 

3.1 Surface Photometry 

The KM models require that the stellar mass spectrum be broken up into a number of 
mass classes with mean stellar mass rrii (see Table 4) each of which is assumed to have an 
energy and angular momentum per unit mass (E and J, respectively) distribution function 
given by 

fi(E, J) OC e -iJ/(^r a )] 2 [e -A t E _ 1]; (1) 

where v s is the scale velocity, W is the reduced gravitational potential, and r a is the 
anisotropy radius beyond which stellar orbits become increasingly radial. For a thorough 
discussion see Da Costa & Freeman (1976) and Gunn & Griffin (1979). To fit the surface 
photometry, Wo determines the shape of the projected luminosity density distribution and 
two parameters, the scale radius r s , and the scale luminosity, are used convert this model 
density to observed units, pc and Lvq pc -2 , respectively. 

The Ai are chosen to be proportional to the mean stellar mass of mass class i, which 
approximates equipartition of energy in the cluster center. Equipartition requires that 
rriiO-f = rrijO-j for all a condition which is not met in the KM formulation due to the 
differing effects of the tidal cut-off on the velocity dispersions of the different mass classes 
(Pryor et al. 1986). The effect mimics incomplete energy equipartition; the high mass 
stars maintain too much kinetic energy relative to the low-mass stars. This is qualitatively 
similar to the Fokker-Planck evolutionary models of Inagaki and Saslaw (1985) in which 
the different mass classes cease to interact after an initial period of energy exchange be- 
cause the lower mass stars are now located at larger radii than the high-mass stars. The 
equipartition/mass segregation dynamics in globular clusters is a complex issue which is 
certainly not fully addressed in the KM formulation. However, by employing both single- 
mass models (no mass segregation, constant mass-to-light ratios) and multi-mass models 
(which probably predict too much mass segregation, see Pryor et al. 1986), we explore the 
range in parameter space within which the true cluster mass segregation state should lie. 

To facilitate the comparison of the KM models, which yield mass density profiles, to 
the SBPs, we need to invoke a stellar mass-luminosity relationship. Bolte (1987) has pub- 
lished a CCD color-magnitude diagram (CMD) of NGC 362. He applied isochrones from 
VandenBerg & Bell (1985), which have heavy element abundances which scale as the solar 
elements, to his photometry. Since that time, Dickens et al. (1991) have carried out high 
resolution echelle spectroscopy of eight member red giants obtaining [Fe/H] = —0.98 ±0.07. 
They also determined that [(C+N+0)/H] = —0.89 ±0.05, slightly overabundant compared 
to solar. We chose to use the stellar models of Bergbusch & VandenBerg (1992) having 
[Fe/H] = -1.03 and [O/Fe] = 0.5. These models have the advantage of extending from the 
tip of the red giant branch (RGB) down to 0.15 M©, but are probably overabundant in 
C,N, and O. We show isochrones produced from this model for ages r = 12, 13 and 14 
Gyr, E(B-V) = 0.0 mag, and (m-M)o = 14.77 mag overlaid on the main-sequence ridge 
line from Bolte (1987) in Fig. 3. We also show one of Bolte's better fits having [M/H] 
= -0.79 and r = 14 Gyr, E(B-V) = 0.0 mag, and (m-M) = 14.60 mag. In the inset 
we have plotted V band luminosity versus stellar mass. The non-existent reddening that 
these isochrones favor is somewhat inconsistent with the value of E(B-V) = 0.04 mag that 
was found by both Burstein & Heiles (1982) from the foreground HI column density and 
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Harris (1982) from UBV photometry. A re-examination of Fig. 3 reveals that it would not 
be possible to reconcile the above isochrones with anything more than a minute amount 
of reddening, although, as Bolte points out, color uncertainties of a few hundredths may 
exist in the photometry calibration and/or the model color-temperature relationship. 

As mentioned, the Bergbusch & VandenBerg isochrones extend to the tip of RGB. 
However, the horizontal branch (HB) stars, despite their short lifetimes, can contribute a 
significant fraction of the total cluster luminosity and even more to the central luminosity. 
The CMD of Harris (1992), demonstrates that the majority of the cluster HB stars are 
red and range in V magnitude from 15.0 to 15.5 with a mean around 15.25 corresponding 
to a mean absolute magnitude of approximately My = 0.5 mag. The main problem with 
attempting to incorporate the HB stars is estimating their lifetimes. Fortunately, the 
HB evolution for stars corresponding to the same model parameters as the Bergbusch & 
VandenBerg isochrones is described in Dorman (1992). A typical lifetime from helium flash 
to core helium exhaustion is 100.0 Myr and is fairly insensitive to the zero-age HB mass. 
By comparing this to the RGB timescales in VandenBerg (1992) we found that 100.0 Myr 
is about the time it takes for a star to evolve up the RGB from the HB level to the RGB 
tip. Therefore, we would expect stars of a similar mass range, approximately 0.001 M© to 
occupy both these evolutionary phases. Using this approach we found that the HB stars 
contributed about 10% of the total cluster luminosity despite contributing less than 0.1% 
of the mass and that the stars in bin 15 of Table 4 contributed 55 - 65% of the luminosity 
and slightly over 1% of the mass. Since post main sequence evolutionary timescales and 
mass luminosity relationships are somewhat uncertain, the large contribution of these stars 
to the cluster luminosity presents a problem in the calculating of synthetic cluster M/L's. 
We will return to this problem later. 

We chose a power-law mass function with a flattening at the faint end: 

(j)( m ) = m" (x+1) dm m> O.3M , (2) 

4>{m) = m dm m < O.3M . (3), 

consistent with Pryor et al. 1989, and 1991. This is similar to what has been seen in 
the galactic disk (Miller & Scalo 1979) but is opposite to what has been derived from 
the luminosity functions of three globular clusters, M13, M71, and NGC 6397, which 
exhibit a steepening at the low-mass end (below 0.4 M ) (Richer et al. 1990). However, 
these luminosity functions were measured at large projected radii where one would expect 
an enhanced number of low-mass stars resulting from cluster dynamical evolution (i.e., 
mass segregation). Cluster environmental effects (i.e. Galactic tidal fields, bulge and disk 
shocking) also tend to have more drastic effects on the outer cluster evolution, although in 
this case one would expect a reduction in low-mass stars. It is fair to say that it is dangerous 
to make conclusions regarding global mass functions based on photometry of stars at large 
projected radii. We will show that it is not possible to reconcile the population and KM 
dynamical M/L's while using a mass function that has a steepening at the low-mass end. 

Remnants were treated in the following manner: stars with initial masses of 0.87 - 1.5 
M , 1.65 - 4.0 M and 4.0 - 8.0 M become white dwarfs with masses of 0.5 M , 0.7 M 
and 1.2 M , respectively (Pryor et al. 1991). These objects are added to the corresponding 
mass bins. More massive stars, which have presumably evolved into neutron stars or black 



9 



holes, are assumed to be ejected from the cluster. This is in agreement with the large 
velocities of many times the cluster escape velocity observed for these objects in the field. 
There is gathering evidence from the many observations of millisecond pulsars (particularly 
in 47 Tuc, Manchester et al. 1991), however, that globular clusters may be more adept 
at holding on to their neutron stars than previously thought. A possible alternative, 
accretion-induced collapse, a method of manufacturing neutron stars through accretion 
on to a white dwarf, has also been suggested (Bailyn & Grindlay 1990). Favoring this 
latter scenario is the difficulty in reconciling the large numbers of millisecond pulsars with 
the small numbers of observed low mass X-ray binaries, the more traditional projenitors 
of millisecond pulsars. Clearly, these considerations would greatly complicate dynamical 
models, particularly for flat mass functions for which the remnant mass classes become 
more significant. As we will show, the addition of a significant number of heavy remnants, 
like a steepening of the low-mass end of the mass function, causes a deterioration in the 
agreement between the dynamical and population M/L's. 

3.2 Radial Velocities 

The mass of a KM model is given by 

M = 9r^l I ^ 2dr (4) 

Illingworth (1976) where v s is the scale velocity. The run of a 2 { (r) and of ^(r) are deter- 
mined from 

I\a t \<W(r) M a ^ W ) a l d3 ^ 

where, o^ = UiCosO or o^sin^ for o r ^ or <Jt,i, respectively, and the i subscript refers to the 
i th mass class. Comparisons were made between the observed velocities and scaled model 
velocity dispersions projected along the line of sight, 

2 ( , 2 p t (r)[(r 2 - R 2 )a 2 jt (r) + R 2 a 2 % (r)}dr 

ap > i[H) H(R)J R r{r 2_ R 2y/2 > W 

(Binney and Tremaine 1987, p. 208), where \ii is the surface density of the i th mass class. 
The optimal velocity scaling and mean velocity were derived using the maximum likelihood 
technique outlined in Gunn and Griffin (1979). 

A serious problem in mass determinations is contamination from binary and non- 
member stars, both of which tend to increase the mass estimate. While non-members 
are a virtual non-issue due to the high systemic radial velocity of NGC 362, binaries can 
be especially problematic with only about one quarter of the stars having multiple radial 
velocity measurements. We attempted to deal with this problem using the technique 
discussed in Fischer et al. (1993). First, using the entire data set, excepting the known 
variables, the optimal v s and v ave are determined using equation 7 for all the KM models 
described below. For every star the parameter 



(r,t),i\ r ) - f TTZ ' 
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was tabulated. The largest Si obtained was approximately 2.80. In order to determine 
the likelihood of having such a value, we performed Monte Carlo simulations of the radial 
velocity data. We started with the known projected radii (Ri) of the program stars. The 
true radius is in the range R < r < r t , where r t is the tidal radius. If x is the displacement 
from the mean cluster position along the line-of-sight such that r = \JR 2 -\- x 2 then the 
probability that the star is at x is 

p(a;)~ P k(VR 2 + x 2 ). (8) 

A three-dimensional position along with corresponding model-dependent radial and tan- 
gential velocities were drawn at random from their respective probability distributions. 
The velocity component along the line-of-sight was then determined, and an error term, 
drawn from a Gaussian distribution with standard deviation equal to the velocity error, as 
tabulated in Table 2, was added. This process was repeated, producing 1000 sets of data, 
each with a given mass and r a and the same projected positions and velocity measurement 
errors as the original data set. Finally the maximum likelihood technique was applied to 
each of the artificial data sets and the maximum 5{ were recorded. We found that a value 
of Smax > 2.80 occurred approximately 70% of the time implying that there is probably 
no high residual velocity contamination. 

3. 3 Results 

The fitted KM parameters for 3 r s < r a < oo and 0.0 < x < 2.0 are shown in Table 
5. Column 1 is the anisotropy radius, column 2 is the mass function slope (the first four 
rows are single-mass models), column 3 is the reduced central potential, column 4 is the 
scale radius, column 5 is the ratio of the tidal radius to the scale radius and columns 6 and 
7 are the reduced chi-squared for the fit {y = 22 degrees of freedom) and the probability 
of exceeding this value, respectively. These probabilities are based on 1000 Monte Carlo 
simulations for each parameter set, each using a surface profile generated from the best fit 
model with errors drawn from the uncertainties shown in Table 1. Hence, the probabilities 
are somewhat dependent on the accuracy of the photometric uncertainties and for this 
reason it is safer to view them in the relative sense. 

The values of v s are shown in column 8 of Table 5. A goodness-of-fit statistic 

^ ("JojU + v err,k) 

was generated for each model and is shown in column 9. The distribution of this statistic 
can be extracted from the radial velocity Monte Carlo simulations described above. We 
find that the ( 2 are distributed around N, the number of radial velocity measurements. 
Values of ( 2 less than N tend to indicate that the model is too flat as a function of 
projected radius for the radial velocity data, while high values imply too steep a model. 
The width of the ( 2 distribution appears to be dependent on the anisotropy radius but 
this is probably due to the higher values of v s for the anisotropic models. These effects are 
demonstrated in Fig 4. The solid histogram in the upper panel is the distribution of £ 2 
when one fits a model with r a = 3r s to data with an isotropic distribution function while 
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the dashed histogram has r a = 3r s for both model and data. Similarly, the solid histogram 
in the lower panel represents an isotropic model and data with r a = 3r s while the dashed 
histogram has isotropic orbits and data. Column 10 of Table 5 shows the probability of 
exceeding the observed |£ 2 — N\ assuming that the cluster velocities are specified by the 
model parameters indicated and have the uncertainties tabulated in Table 2. The greater 
this probability the higher the likelihood that the cluster velocities are drawn from the 
specified distribution. 

Derived parameters for the KM model are displayed in Table 6: columns 1 and 2 
specify the anisotropy radius and mass function slope, while columns 3 and 4 contain the 
central luminosity density and total cluster luminosity, respectively. Columns 5 and 6 are 
the central mass density and total cluster mass, respectively. Columns 7 and 8 are the 
central and global population M/L's given by 



where l(m) is the luminosity of a star of mass m given by the Bergbusch & VandenBerg 
stellar models described above, and mi and m u are the lower and upper mass cut-offs, 
respectively. 

The corresponding dynamical M/L's are in columns 9 and 10, respectively. Monte- 
Carlo simulations of the radial velocity data, as described above, were use to determine the 
uncertainties implicit in the maximum likelihood technique and to search for any possible 
systematic effects. We noticed that the maximum likelihood method resulted in scale 
velocities that were biased systematically too low but only by a few percent so they have 
not been corrected. This bias tends to be larger when the number of radial velocities is 
smaller. 

First a comment about the single-mass models which assume that stars of all masses 
have the same radial distribution. We obtain values of c and r s consistent with both 
Illingworth & Illingworth (1976) and Tucholke (1992) (both used only isotropic models) 
indicating that our luminosity profile is very similar to their starcounts. Assuming isotropy, 
we find a higher central velocity dispersion (although consistent within uncertainties) than 
the value of 7.8 km s _1 which was obtained by Illingworth (1976) from integrated spec- 
troscopy of the cluster center (this value was corrected from the observed velocity dispersion 
of 7.5 km s _1 assuming an isotropic model). Correspondingly, our isotropic single-mass 
mass estimate is about 30% higher which puts us right at the upper limit of Illingworth's 
la uncertainty. Our luminosity estimate is 17% lower, resulting in a M/L = 1.40 ± 0.2 
M©/Ly ©which is 55% higher than his, a difference significant at the 1.4cr level. The clus- 
ter has a mean systemic velocity, independent of model parameters, of 223.5 ± 2.0 km 
s _1 where the uncertainty reflects the accuracy of our zero-point and includes the actual 
velocity errors in the standard system itself. This is in excellent agreement with the in- 
tegrated radial velocity of Dubath et al. (1993) who obtained 223.20 ± 0.2 km s _1 , and 
may indicate that our velocity zero-point is more secure than 2.0 km s _1 . Fig. 5 shows 
the surface brightness data from Table 1 along with the isotropic and r a = 3r s single-mass 
models. Also shown as a dotted line is a typical stellar profile having FWHM = 1.80" 
which is about 14% of the cluster core radius. Therefore, the core is easily resolved and 
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the seeing will have negligible impact on its measurement. Finally, we have plotted a 
power-law with slope equal to -1 as a long-dashed line. NGC 362 has been classified as 
a possible post core-collapse cluster (PCC? in their designation scheme) by Chernoff & 
Djorgovski (1989). The criteria that these authors used to determine if a cluster is PCC 
are: a small core, and "... an unambiguous, extended, slope ~ —1 power-law section of 
the surface brightness profile near the center...". Neither of these criteria are applicable to 
our data and we conclude that NGC 362 is not a post core-collapse cluster. Unfortunately 
Chernoff & Djorgovski did not include a plot of their SBP of NGC 362. 

For the multi-mass models, the trends in the goodness-of-fit to the SBP with r a 
and x can be categorized as follows: 1) High quality fits can be found for the whole x 
and most of the r a range shown in Tables 5 and 6, but not for all combinations of the 
parameters. 2) For low values of x the best fits occur for moderate to high amounts of 
anisotropy. 3) For high values of x only the isotropic models provide good fits. This 
is easily explained; as the mass function steepens, one adds more low-mass stars to the 
cluster, which, because of equipartition of energy migrate to large radius giving the cluster 
a more extended appearance. This mimics the effect of decreasing r a for shallower mass 
functions which also gives the cluster a more extended appearance due to radial orbits 
traveling to larger cluster-centric radii than equally energetic circular orbits. Overall, the 
SBP is more consistent with shallow mass functions and moderately low levels of anisotropy 
for the range of parameter space sampled here. 

The trends are much simpler for the radial velocity modeling in that highly anisotropic 
orbits are always favored, regardless of mass function. In fact, it is possible to rule out all 
but the x = 0.0 isotropic models at a very high level of confidence. The best overall fits 
to both the kinematics and the SBP occur for x = 0.0 — 0.5 and r a = 5 — 10 It should be 
pointed out, however, that while the x 2 result from the fitting of the SBP, the ( 2 arise from 
the fitting of radial velocities to velocity models whose shape is fixed by the SBP leaving 
only a velocity scaling and mean velocity to be determined. In other words, we have seen 
that the radial velocities are most consistent with anisotropic distribution functions which 
have radial velocity dispersion profiles which fall off more rapidly at small radii than their 
isotropic counterparts. A similar effect can be accomplished by decreasing the value of r s 
from the one most consistent with the SBP. Therefore, what we are really showing here 
is not that the radial velocity data are only consistent with anisotropic models, but that 
the only way to get consistency for all the data is to invoke moderately anisotropic models 
and shallow mass functions. 

One apparent problem with this set of KM models is that the central and global 
population M/L's are all substantially higher than the dynamic M/L's, and the situation 
is worse for the best-fitting models with x = 0.0 — 0.5 and r a = 5 — 10. These models have 
both global and central population M/L's which are 1.5 - 2.5 times larger than the dynamic 
estimates. How much of a problem is this? The population M/L's are very sensitive to how 
the RGB and HB are populated and different assumptions can lead to very different results 
as can be seen in Pryor et al. (1991) who obtain global M/L's that are 25 - 50% lower 
than our tabulated values. Such a reduction would provide significantly better agreement 
between the population and dynamical M/L's, especially for x = 0.5. Therefore, we have 
arbitrarily increased the number of stars in bin 15 by a factor of two. This is equivalent 
to doubling the (uncertain) lifetimes on the RGB and the HB. Aside from the uncertainty 
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in these lifetimes, another justification is that the population models do not include all 
the possible evolved stars (i.e. asymptotic branch stars) or the possible existence of blue 
stragglers many of which may be hiding in the core as is the case for 47 Tuc (Parasce et al. 
1991). The results are shown in Tables 7 and 8 which have the same format as Tables 5 and 
6, respectively. There is very little change in most of the derived and fitted parameters 
or the quality of the fits for a given set of parameters. However, one can see that the 
desired effect has been achieved, the values of the global and central population M/L's 
are reduced and the agreement, particularly for x = 0.5 — 1.0, is considerably improved. 
Clearly, with similar manipulations one could achieve good agreement with virtually any x 
and this tends to make comparisons between dynamic and population M/L's a fairly weak 
constraint. The situation could be greatly improved with complete number counts of the 
evolved stars; a difficult task requiring high resolution imaging. 

To summarize, the best fits to both the kinematic data and surface photometry occur 
for a flat mass function with moderate amounts of anisotropy. The cluster luminosity is 
around 1.70 ±0.1 x 10 5 L VQ and the mass ranges from approximately 2.5 — 3.5 x 10 5 M 
for a global mass-to-light in the range M/L v = 1.5 - 2.0 M /Ly with a mass function 
near x = 0.5. Unlike the findings of Richer et al. (1990) for 3 other clusters, our data 
is inconsistent with a sharp upturn in the mass function at the low-mass end. Such an 
upturn would result in a large increase in the population M/L requiring a large increase 
in the number of giants or a decrease in the number of stellar remnants both of which 
would be hard to justify. It would also be possible to get good agreement between the 
SBP and the kinematics for mass functions less than zero and isotropic orbits (not shown 
in tables) However, this results in a large divergence between population and dynamic 
M/L's. Including large numbers of heavy remnants (black holes, neutron stars) results in 
a similar problem. 

Recently Capaccioli et al. (1991) have shown that there exists a fairly strong correla- 
tion between the slope of cluster mass functions and both the clusters' distance from the 
Galactic center (Rgc) and above the Galactic disk (Z). There are several reasons why one 
would expect shallower cluster mass function near the Galactic disk and center. Perhaps 
the most compelling is heating of the cluster through disk and tidal shocking, although 
Galactic tidal fields will also play a role. The effect is to preferentially evaporate low-mass 
stars, which tend to be more loosely bound to the cluster, resulting in a present day mass 
function (PDMF) which is flatter than the initial mass function (IMF). NGC 362 with 
(m-M)o = 14.77 mag, / = 301.53°and b = -46.25°, has R GC = 9.6 kpc and Z = 6.5 kpc 
(assumes the distance from the sun to the galactic center is 8 kpc). Therefore, NGC 362 
has a similar value of x and Rgc to M92 and hence agrees fairly well with the weaker 
x - Rgc correlation. However, it tends to be an outlier in the stronger x - Z correlation 
(assuming it was extrapolated to 6.5 kpc) requiring a value of about x = 1.0 for consis- 
tency. NGC 362 has a larger value of Z than all but one cluster in the Capaccioli et al. 
compilation; the largest is M3 with Z = 10.2 which also lies below the correlation. The 
evidence from these two high Z clusters may indicate that the relationship is flattening. 
Perhaps this implies that NGC 362 is sufficiently far from the disk that tidal shocking is 
relatively unimportant to its dynamical evolution and that its mass function has evolved 
to a lesser degree than those of the low Z clusters. A more interesting correlation to search 
for, in terms of disk and bulge shocking, would be the one between x and various cluster 
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orbital parameters, particularly the clusters' perigalacticon. The mass functions should be 
more sensitive to these parameters, as they more truly reflect the environmentally driven 
aspects of cluster evolution (i.e. see Aguilar et al. 1988). 

3.4 Relaxation Timescales 

As mentioned above, a real cluster cannot achieve complete equipartition of energy 
regardless of its age. However, since, through the use of KM models, we are assuming a 
large degree of mass segregation, it is necessary to justify this in terms of the cluster's 
evolutionary state. Two relevant timescales are the oft-quoted central relaxation time 

t ro = (1.55 x 10V) (j^pr) (^) Pog(0-5M/ (m))]" 1 = (0.05 - 1.0) x 10 8 yr 

(11), 

(Lightman and Shapiro 1978) and the half mass relaxation time 

t rh = (8.92x10V) (jT^y (H) Pog(0.4M/ (m))]" 1 = (0.7-2.3) x 10 9 yr 

(Spitzer and Hart 1971). 

The ranges in relaxation time correspond to the parameter sets described above. It 
appears that the system is largely relaxed therefore necessitating the use of multi-mass 
models. 

4. Conclusions 

In this paper we have examined the internal dynamics of the globular cluster NGC 
362 using V band CCD images and echelle spectra of 208 member red giants. 

1) The new oxygen-enhanced isochrones of Bergbusch & VandenBerg (1992) with [Fe/H] 
= -1.03 were applied to the fiducial main sequence and sub-giant branch (Bolte 1987) 
of NGC 362 yielding an age of 13 Gyr, an absolute distance modulus (m-M)o = 14.77 
mag, and reddening E(B-V) = 0.0 mag. 

2) A CCD surface brightness profile (SBP) extending out to approximately 20 pc was 
constructed from a mosaic of four V band images. After determining that NGC 362 
is not post core-collapse, single- and multi-mass King Michie models were applied 
to the data for a range of mass function slopes 0.0 < x < 2.0 and anisotropy radii 
3.0 < r a (r s ) < 00. The quality of the fits varied as a function of both these parameters; 
for small values of x, anisotropic models were favored while for steeper mass functions 
isotropic orbits provided superior fits. The derived total cluster luminosity is 1.70 
±0.1 x 10 5 L VQ . 

3) A total of 285 stellar spectra of 215 stars were obtained; 50 were observed twice, one 
observed three time, one observed ten times, and one observed eleven times. Three 
stars were obvious non-members and four stars had strong evidence for radial velocity 
variability which we conclude is likely due to the presence of binary stars. The implied 
cluster binary fraction is 0.15 for circular orbits or 0.27 for an orbital distribution 
function /(e) = e. This is a considerably higher binary fraction than previous cluster 
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studies have found. While it may indicate that NGC 362 has a high binary fraction, 
it might also be a result of the way we selected candidates for repeat measurements. 
Furthermore, this fraction is based on only two measurements of each of the candidate 
binaries. 

4) The 208 remaining stars with radial velocity measurements were incompatible with 
King-Michie models having isotropic velocity dispersion and luminosity profiles con- 
sistent with the SBP. Therefore, the best agreement with both the kinematic data and 
the SBP were for shallow mass functions x = 0.0 — 0.5 and intermediate amounts of ve- 
locity anisotropy. In this best-fit range, the derived cluster mass was M = 2.5— 3.5x 10 5 
M for a global mass-to-light ratio of M/Ly = 1.5 - 2.0 M q /Lvq- 

5) The values of x and Z, the cluster height above the Galactic disk, do not lie on 
the correlation (assuming it was extrapolated to Z = 6.5 kpc) between these two 
parameters found by Capaccioli et al. (1991). NGC 362 has a mass function which is 
too flat, possibly implying that the cluster evolutionary drivers (i.e. disk and bulge 
shocking, tidal fields) play a relatively small role at its location in the Galaxy. 

6) Our value of x differs from the sharp turn-ups in mass functions derived from some 
deep luminosity function work reported in Richer et al. (1990) for the clusters M13, 
M71, and NGC 6397. 

7) In the future we will attempt non-parametric modeling of the cluster radial velocities 
as described by Merritt (1992). 
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Fig. 1- Mean radial velocity vs. projected radius (upper panel) and versus position angle 
(lower panel) for the 212 member stars. The solid lines are the mean velocity, v = 223.5±2.0 
km s _1 . 

Fig. 2- The difference in median velocity for stars on either side of an axis at the specified 
position angle. Also shown is the best fit sine function with the parameters shown at the 
top of the plot. 

Fig. 3- A plot of the fiducial main sequence ridge (dotted line) from Bolte (1987). The solid 
lines are isochrones from Bergbusch & VandenBerg (1992) for [Fe/H] = -1.03, Y=0.2368, 
(m-M) = 14.77 mag, and E(B-V) = mag with ages of r = 12, 13, and 14 Gyr. The 
long-dashed line is an isochrone from VandenBerg & Bell (1985) with for [Fe/H] = -0.79, 
Y=0.3, (m-M) = 14.60, and E(B-V) = mag with an age of r = 14 Gyr. 
Fig. 4- The solid histogram in the upper panel is the distribution of £ 2 when one fits 
a model with r a = 3r s to data with an isotropic distribution function while the dashed 
histogram has r a = 3r s for both model and data. Similarly, the solid histogram in lower 
panel represents and isotropic model and data with r a = 3r s while the dashed histogram 
has isotropic orbits and data. 

Fig. 5- The surface brightness data (squares) along with the isotropic (solid line) and 
r a = 3r s (dashed line) single-mass models. The dotted line is a typical stellar profile and 
the long-dashed line is a power law relationship with slop equal to -1. 
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0.77 


0.01 


1 


H2403 01 


•D3 


r il 


21 


i vj. <jvj. 


97 


7 


1 83 8^ 

lOO. OtJ 


£1 8 

O J- . o 


47871 

11 OI -1 


■ "JvJ^I 


226. 


4fi 


+ 


1 


Q3 

. .JO 












H2404 01 




. ij^ 


46 


i VJ. o\j 


22 


u 


9DD 43 


si n 


47871 

11 OI -1 , 


c ;fi44 

■ vjvjll 


93^ 

£jOO 


OCf 


+ 


1 


^6 

. OK) 
































48605 


.6271 


235. 


.15 


± 


2 


.26 


235.45 


± 


1.28 


0.03 


1 


H2410 01 


:03: 


50 


.50 


-70:49: 


22, 


4 


200 65 


62.9 


47871. 


5666 


221 


23 


± 


2 


.18 












H2411 01 


:03: 


40 


.21 


-70:49: 


:41. 


.8 


146.82 


60.6 


47871. 


.5686 


223. 


.98 


± 


1 


.65 












H2417 01 


:03: 


32 


.82 


-70:48: 


:53. 


,1 


151.59 


37.1 


47871. 


.5742 


222. 


15 




1 


.90 












H2418 01 


:03: 


32 


.13 


-70:48: 


40, 


2 


160.17 


33.4 


47871. 


.5780 


218. 


.37 




1 


.53 












H2419 01 


:03: 


36 


.50 


-70:48: 


:47. 


.6 


167.29 


40.9 


47871. 


.5761 


232 


.23 


± 


1 


.67 
































48605 


.6354 


233. 


.06 


± 


1 


.35 


232.73 




1.05 


0.15 


1 


H2423 01 


:03: 


34 


.89 


-70:48: 


07, 


,7 


194.91 


31.5 


47871. 


.5799 


223. 


10 


± 


o 


.92 












H2431 01 


:03: 


24 


.04 


-70:47: 


50, 


4 


189.80 


14.7 


47871. 


.5813 


224. 


.74 


± 


o 


.77 












cl 


01 


•03- 


1 3 

± o 


R4 


-70^1 


1 3 
±o . 


7 


1 Q Q7 


1 88 

lOO. JU 


47867 

1 / OVJ / 


^987 


220. 


. vJU 


± 


n 


71 

. 1 -i 
































47868 


.5365 


219. 


70 


± 


o 


.70 


220.12 


± 


0.50 


0.74 


1 


aa 


01 


:03: 


15 


.46 


-70:51: 


06, 


5 


13.91 


154.44 


47868. 


.5313 


204. 


.70 


± 


1 


.52 












ab 


01 


:03: 


15 


.47 


-70:51: 


09, 


.6 


16.78 


158.86 


47868. 


.5328 


216 


.94 


± 


1 


.04 












tic 


01 


:03: 


14 


.67 


-70:51: 


15 


o 


21.16 


174.26 


47868. 


.5346 


227. 


01 


± 


o 


.81 












ad 


01 


:03: 


14 


.22 


-70:51: 


15 


,8 


21.85 


180.26 


47868. 


.5383 


233 


12 


± 


o 


.82 












ae 


01 


•03- 


1 3 

J- O 




-TO- "SI 

1 VJ. rJ -1 


22, 


o 


98 71 


1 8Q 07 


47868 


^409 


218. 


^9 


± 


n 


7D 

. 1 VJ 












af 


01 


:03: 


15 


.13 


-70:51: 


26, 


.8 


33.14 


172.41 


47868. 
47874 


.5420 
.5410 


208. 
208. 


11 

.79 


-H -H 


o 
o 


.74 
.91 
































48607 


.5508 


208. 


.92 


± 


o 


.82 


208 56 


± 


0.47 


0.63 


2 


a g 


U -L 


:03: 


15 


.01 


-70:51: 


32 


5 


38.74 


174.39 


47868. 


.5437 


223. 


.09 


± 


1 


.04 
































47874 


^497 


222. 


D8 

. uo 


± 


1 


97 


999 R8 


+ 


n 8D 

VJ. OVJ 


38 

VJ. oo 


1 


ah 


01 


:03: 


16 


.05 


-70:51: 


22, 


.6 


30.00 


162.74 


47868. 


.5455 


217, 


.46 


± 


1 


,23 
































47874 


.6135 


215 


75 


± 


1 


.33 


216.67 


-t 


0.90 


0.89 


1 


aj 


01 


:03: 


16 


.76 


-70:51: 


11 


.6 


21.57 


144.92 


47868. 


.5474 


232. 


,15 


± 


1 


.15 












ak 


01 


:03: 


11 


.06 


-70:51: 


02, 


1 


17.64 


242.48 


47868 


.5541 


227, 


15 


± 


1 


.12 












am 


01 


•03- 

.VJO. 


i n 


3D 


i vj. <j _i . 






22.99 


937 4Q 


478R8 

1 / OVJO 


■ OOOO 


228. 


8fi 
ovj 


± 


1 


95 












an 


01 


•03- 


11 


8fi 
. oij 


1 VJ. rJ .1 


•09 


u 


14 ^fi 


933 55 

iOO . oo 


47868 

1 I OvJO 


■ OOOO 


223. 


3^ 

OO 


± 


n 

VJ 


Q7 












ap 


01 


:03: 


12 


.63 


-70:50: 


50, 


o 


8.85 


296.49 


47868. 


.5633 


223. 


59 


± 


1 


.43 












aq 


01 


:03: 


13 


.46 


-70:50: 


50, 


1 


5.44 


315.08 


47868. 


.5651 


217 


.22 


± 


1 


.26 












ar 


01 


:03: 


13 


.48 


-70:50: 


52, 


.9 


3.88 


285 68 


47868. 


.5668 


218. 


.86 


± 


1 


.42 












as 


01 


:03: 


.14 


.08 


-70:50: 


:48 


3 


5.70 


352.07 


47868. 


.5688 


224. 


.34 


± 


1 


.08 












at 


01 


•rn 

.uo. 


12 


^6 

. <jvj 


( VJ. .JVJ. 


•37 
.01. 


q 




339 74 


478R8 

1 / OvJO 


^704 

. O i VJ1 


221. 


Q4 


± 


n 


69 
































4 / o / o 




999 
ZZZ. 


no 

.yo 


± 


n 
U 


89 


000 ^0 


1 


n aq 


1 no 


i 
i 


au 


01 


:03: 


12 


.95 


-70:50: 


:35, 


.3 


19.70 


341.20 


47868 


.5721 


208. 


.60 


± 





.92 












av 


01 


:03: 


14 


.94 


-70:51: 


:00. 


.1 


7.05 


150.75 


47868. 


.5743 


229. 


.62 


± 





.68 












aw 


01 


:03: 


;15 


.93 


-70:50: 


:54, 


.0 


8.32 


90.3 


47868. 


.5762 


219 


.67 


± 


1 


.01 












ax 


01 


:03: 


19 


.16 


-70:50: 


:59. 


.5 


24.84 


102.92 


47868. 


.5781 


228. 


20 


± 


1 


.05 












ay 


01 


:03: 


19 


.00 


-70:50: 


:51. 


.2 


23.59 


83.3 


47868. 


.5797 


222. 


15 


± 





.71 












az 


01 


:03: 


17 


.55 


-70:50: 


:41. 


.8 


20.32 


53.2 


47868. 


.5816 


224. 


.64 


± 


1 


.32 












b 


01 


:03: 


13 


.77 


-70:51: 


:09. 


.1 


15.33 


188.68 


47867. 


.5307 


208. 


.12 


± 





.61 












ba 


01 


:03: 


19 


.32 


-70:50: 


:34. 


.0 


31.99 


51.4 


47869. 


.5228 


233 


.67 


± 





.62 












bb 


01 


:03: 


20 


.62 


-70:50: 


:29. 


.8 


39.61 


52.4 


47869. 


.5287 


238. 


13 


± 





.71 
































48607 


.5479 


239. 


,66 


± 





.99 


238.65 


± 


0.58 


1.58 


1 


be 


01 


:03: 


16 


.44 


-70:50: 


:26 


.9 


29.14 


21.8 


47869. 


.5304 


221. 


.21 


± 


1 


.53 













Table 2 (cont.) 
Radial Velocities 
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Star a 5 R PA HJD v r v r \ v Comments 

(2000.0) (2000.0) (arcsec) (°) -2400000 (kms" 1 ) (km s _1 ) 



bd 


01:03:14.99 


-70:50:26.8 


27.40 


7.74 


47869.5321 


233.03 


± 


1. 


.12 












be 


01:03:14.77 


-70: 


:50:14.5 


39.54 


3.78 


47869.5339 


220.51 


± 


1. 


.27 












bf 


01:03:12.90 


-70: 


:50:07.7 


46.72 


351.88 


47869.5357 


224.95 


± 


0. 


,97 












bg 


01:03:11.54 


-70: 


:50:22.7 


33.96 


336.96 


47869.5375 


224.82 


± 


0. 


.91 












bh 


01 


:03: 


:10.24 


-70: 


:50:15.5 


43.20 


332.88 


47869.5392 


226.67 


± 


0. 


.79 












bi 

J 


01 


:03: 


:09.03 


-70: 


:50:20.0 


42.55 


322.92 


47869.5409 


217.02 


± 


0. 


,70 












bk 


01 


:03: 


:09.99 


-70: 


:50:21.8 


38.36 


326.94 


47869.5426 


226.17 


± 


1 


.20 












bm 


01 


:03: 


:05.55 


-70: 


:50:30.9 


48.59 


298.30 


47869.5461 


225.73 


± 


0. 


,84 












bn 


01 


:03 


:05.43 


-70: 


:50:28.2 


50.43 


300.69 


47869.5478 


224.88 


± 


1 


,19 












bp 


01 


:03: 


:06.05 


-70: 


:50:36.6 


43.88 


293.27 


47869.5540 


217.38 


± 


1. 


,67 












bq 


01 


:03: 


:06.69 


-70: 


:50:35.5 


41.49 


296.39 


47869.5557 


225.86 


± 


1. 


,77 












br 


01 


:03: 


:07.22 


-70: 


:50:48.5 


34.97 


278.95 


47869.5575 


225.00 


± 


1 


,33 












bs 


01 


:03: 


:09.88 


-70: 


:50:48.2 


22.21 


284.99 


47869.5592 


224.57 


± 


2. 


,15 












bt 


01 


:03: 


:05.65 


-70: 


:50:59.0 


42.57 


263.17 


47869.5627 


223.41 


± 


0. 


.71 












bu 


01 


:03: 


:04.78 


-70: 


:50:58.7 


46.79 


264.16 


47869.5643 


228.49 


± 


0. 


,87 












bv 


01 


:03: 


:06.06 


-70: 


:51:09.2 


43.04 


249.23 


47869.5661 


222.72 


± 


0. 


,96 












bw 


01 


:03: 


:09.62 


-70: 


:51:00.4 


23.63 


254.15 


47869.5680 


222.17 


± 


2. 


,97 












bx 


01 


:03 


:11.69 


-70: 


: 50:51.1 


12.87 


282.79 


47869.5698 


222.19 


± 


1 


,35 












by 


01 


:03: 


:09.15 


-70: 


:51:16.7 


33.83 


227.74 


47869.5715 


236.62 


± 





,60 




























47873.5612 


234.79 


± 





,80 


235.96 


± 


0.48 


3.35 


1 


bz 


01 


:03: 


11.65 


-70: 


:51:25.9 


34.40 


201.74 


47869.5733 


228.70 


± 


1. 


93 












c 


01 


:03 


:13.00 


-70: 


:51:14.3 


21.25 


196.69 


47867.5325 


209.72 


± 


1 


,14 












ca 


01 


:03: 


:11.57 


-70: 


:51:37.1 


45.11 


196.93 


47869.5750 


231.81 


± 


0. 


,91 












cb 


01 


:03: 


10.19 


-70: 


:51:35.1 


45.72 


205.83 


47869.5767 


216.34 


± 


2. 


,55 












cc 


01 


:03: 


:12.15 


-70: 


:51:41.3 


48.45 


192.25 


47869.5785 


221.65 


± 


0. 


,80 












cd 


01 


:03: 


13.26 


-70: 


:51:33.4 


39.74 


186.97 


47869.5804 


222.85 


± 


1. 


,23 












ce 


01 


:03: 


:07.80 


-70: 


: 50:11. 3 


53.14 


323.36 


47870.5276 


220.39 


± 


1 


23 












cf 


01 


:03: 


:07.37 


-70: 


: 50:11.8 


54.04 


321.25 


47870.5293 


212.28 


± 


0. 


,96 




























48607.5436 


212.15 


± 


1 


.10 


212.22 


± 


0.72 


0.01 


1 


ce 

o 


01 


:03: 


:09.07 


-70: 


:50:53.4 


25.45 


271.23 


47870.5310 


218.64 


± 


0. 


.97 












ch 


01 


:03 


:12.24 


-70: 


:51:33.5 


40.76 


193.97 


47870.5347 


221.22 


± 


1 


.16 












cj 


01 


:03 


:17.56 


-70: 


:51:18.4 


29.40 


146.26 


47870.5365 


220.90 


± 


1 


,98 












ck 


01 


:03 


:16.53 


-70: 


:51:02.1 


13.91 


125.88 


47870.5399 


223.56 


± 


1 


.53 












cm 


01 


:03: 


:18.84 


-70: 


:51:14.8 


30.77 


132.66 


47870.5418 


224.94 


± 


0. 


,97 




























47874.6171 


223.80 


± 


1. 


,06 


224.42 


± 


0.72 


0.63 


1 


cn 


01 


:03: 


:18.91 


-70: 


:51:26.1 


39.52 


144.46 


47870.5436 


231.95 


± 


1 


,26 




























47870.5493 


232.53 


± 


1. 


,10 


232.28 


± 


0.83 


0.12 


1 


cp 


01 


:03: 


: 17.47 


-70: 


:51:30.0 


39.40 


156.22 


47870.5455 


221.31 


± 


1, 


,91 












cq 


01 


:03: 


:15.96 


-70: 


:51:36.3 


43.19 


168.70 


47870.5474 


227.15 


± 


1 


,63 












cr 


01 


:03: 


:20.83 


-70: 


:51:22.0 


42.87 


130.88 


47870.5513 


220.32 


± 


0. 


,89 




























47874.6153 


220.43 


± 


1. 


.07 


220.36 


± 


0.68 


0.01 


1 


cs 


01 


:03 


:21.35 


-70: 


:51:20.7 


44.04 


127.42 


47870.5577 


219.62 


± 


1, 


.21 












ct 


01 


:03: 


:21.80 


-70: 


:51:27.0 


49.76 


131.64 


47870.5599 


219.06 


± 


1 


30 












cu 


01:03:23.55 


-70: 


:51:21.0 


53.20 


120.58 


47870.5618 


205.17 


± 


1, 


.16 




























48605.6023 


206.70 


± 


1. 


.57 


205.71 


± 


0.93 


0.61 


1 


cv 


01:03:23.65 


-70: 


:51:08.5 


48.53 


107.46 


47870.5637 


223.82 


± 


0. 


,74 












cw 


01:03:24.61 


-70: 


:51:08.9 


53.17 


106.35 


47870.5656 


222.75 


± 


1, 


,37 












cx 


01:03:24.64 


-70: 


:51:01.1 


51.67 


97.9 


47870.5676 


220.97 


± 


2. 


,75 













Table 2 (cont.) 
Radial Velocities 
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Star a 5 R PA HJD v r v r \ v Comments 

(2000.0) (2000.0) (arcsec) (°) -2400000 (kms" 1 ) (km s _1 ) 



cy 


01 


:03: 


:23. 


92 


-70: 


:50: 


:45. 


.1 


48, 


,45 


79.5 


47870, 


.5694 


219, 


.24 


± 


1. 


43 












cz 


01 


:03: 


:21. 


33 


-70: 


:50: 


:49. 


.7 


35. 


,15 


83.0 


47870, 


,5715 


216, 


.59 


± 


2, 


,48 












d 


01 


:03: 


:12. 


.61 


-70: 


:51: 


:02. 


.9 


12. 


,02 


221.86 


47867, 


.5349 


214, 


,11 


± 


0. 


,79 












da 


01 


:03: 


:18. 


.70 


-70: 


:50: 


:32, 


.8 


30. 


,48 


46.0 


47870, 


.5735 


221. 


.44 


± 


1 


,27 












db 


01 


:03: 


:17. 


.34 


-70: 


:50: 


:26. 


,3 


31. 


,58 


28.9 


47870, 


.5752 


216. 


.13 


± 


1. 


,52 












dc 


01 


:03: 


:23. 


91 


-70: 


:50: 


:24, 


.5 


55. 


.97 


58.2 


47870, 


,5770 


221, 


.63 


± 


1. 


,02 












de 


01 


:03: 


:24. 


.31 


-70: 


:50: 


:19 


.2 


60 


.53 


54.9 


47870, 


,5787 


221. 


.46 


± 


1. 


.42 




































47871, 


,5485 


223 


.78 


± 


1, 


.13 


222.88 


± 


0.88 


1.63 


1 


df 


01 


:03: 


:21. 


.70 


-70: 


:50: 


A3. 


.3 


38, 


,23 


73.8 


47870, 


,5808 


226. 


.60 


± 


2, 


,31 












dg 


01 


:03 


:16 


.38 


-70: 


:50: 


:04 


.7 


50, 


.36 


12.0 


47870, 


5831 


224. 


.72 


± 


1 


.25 












dh 


01 


:03: 


:06 


.46 


-70: 


:51: 


AO. 


.0 


59, 


,88 


219.72 


47870, 


,5849 


225. 


,04 


± 


1 


,07 












di 


01 


:03 


:03 


.89 


-70: 


:51: 


:31. 


.9 


63 


,51 


233.28 


47870, 


,5867 


220. 


.39 


± 





,68 












dk 


01 


:03: 


:25. 


,78 


-70: 


:49: 


:47. 


,9 


87. 


.12 


40.7 


47871, 


,5289 


207, 


.16 


± 


2. 


,03 




































48605, 


,6207 


212 


.71 


± 


1. 


,93 


210.08 


± 


1.40 


3.93 


1 


dm 


01 


:03: 


:25. 


.19 


-70: 


:49: 


:58. 


.9 


77. 


.05 


44.4 


47871. 


.5310 


224, 


.78 


± 


0, 


,84 












dn 


01 


:03: 


:23. 


.13 


-70: 


:49: 


:54. 


.6 


73. 


,74 


36.4 


47871. 


,5328 


223. 


.72 


± 


2. 


,03 












dp 


01 


:03: 


:20. 


.01 


-70: 


:49: 


:55. 


,2 


65. 


.26 


25.8 


47871. 


,5344 


230 


.01 


± 


1. 


,47 












dq 


01 


:03: 


:20. 


.54 


-70: 


:49: 


:36. 


.8 


83, 


15 


21.9 


47871. 


,5364 


217. 


.51 


± 


2 


.19 












dr 


01 


:03 


:21 


.97 


-70: 


:49: 


:58. 


,2 


67, 


,50 


34.3 


47871. 


,5389 


223 


.41 


± 


1 


.98 












ds 


01 


:03: 


:25. 


.00 


-70: 


:50: 


:08. 


.5 


69, 


,79 


49.3 


47871, 


5411 


220. 


.58 


± 


2, 


,07 












dt 


01 


:03 


:28 


.37 


-70: 


:50: 


AO. 


,8 


70, 


,77 


79.3 


47871. 


,5549 


222. 


.70 


± 


1 


32 












du 


01 


:03: 


:11 


.09 


-70: 


:49: 


A3. 


,4 


72, 


23 


347.60 


47872, 


,5686 


221. 


.13 


± 


1 


,26 












dv 


01 


:03: 


:07. 


.76 


-70: 


:49: 


:46 


.2 


74. 


,89 


334.77 


47872, 


.5704 


234, 


.26 


± 


0. 


,69 




































48605, 


,6518 


236 


.18 


± 


1. 


,02 


234.86 


± 


0.57 


2.43 


1 


dw 


01 


:03: 


:03. 


.15 


-70: 


:50: 


:30. 


.2 


59, 


,53 


293.49 


47872, 


.5722 


227, 


.12 


± 


0. 


,71 












dx 


01 


:03: 


:02. 


,78 


-70: 


:50: 


16, 


.9 


67. 


,49 


303.27 


47872, 


.5739 


216. 


.03 


± 


0, 


93 












dy 

J 


01 


:03: 


:01. 


,43 


-70: 


111 


:59. 


.8 


83. 


.12 


310.63 


47872, 


.5756 


220. 


.18 


± 


0. 


,95 




































47873, 


,5629 


220 


.28 


± 


0. 


,62 


220.25 


± 


0.52 


0.01 


1 


dz 


01 


:03: 


:01. 


.12 


-70: 


:50: 


:19 


.3 


73 


.29 


298.19 


47872, 


5773 


212. 


.94 


± 


2, 


.01 




































48605, 


,6535 


213 


.16 


± 


1 


,54 


213.08 


± 


1.22 


0.01 


1 


e 


01 


:03: 


:10. 


.62 


-70: 


:50: 


:53. 


.8 


17, 


.81 


270.47 


47867, 


5367 


233 


.57 


± 


0. 


,72 




































47873, 


,6263 


233 


.17 


± 





,80 


233.39 


± 


0.54 


0.14 


1 


ea 


01 


:03: 


:01. 


,92 


-70: 


:51: 


:17 


.5 


65, 


03 


248.74 


47873, 


,5646 


235. 


.29 


± 


1. 


,68 




































48605, 


,6585 


238 


.98 


± 


1. 


,48 


237.37 


± 


1.11 


2.72 


1 


eb 


01 


:03: 


:00. 


13 


-70: 


:51: 


:00. 


.2 


69. 


.71 


264.83 


47873, 


.5669 


217. 


.24 


± 


1. 


,62 












ec 


01 


:02: 


:57. 


93 


-70: 


:50: 


AO. 


.4 


81. 


,41 


279.55 


47873, 


.5694 


225. 


.53 


± 


1. 


,64 












ed 


01 


:03: 


:08. 


,91 


-70: 


:51: 


:49. 


.6 


61. 


,52 


205.22 


47873, 


,6182 


226. 


.04 


± 


0, 


.89 












edd 


01 


:03: 


:15. 


41 


-70: 


:51: 


:42. 


.6 


48, 


99 


173.26 


47874, 


,5444 


223. 


.69 


± 


1. 


,01 












ef 


01 


:03 


:10 


.94 


-70: 


:51: 


:54, 


,2 


62, 


.40 


195.07 


47873. 


6197 


215. 


.07 


± 





,93 




































48607, 


,5690 


215 


.01 


± 


0. 


,78 


215.03 


± 


0.60 


0.00 


1 


eff 


01 


:03 


:14 


.26 


-70: 


:51: 


:58, 


.0 


64, 


.05 


179.91 


47874. 


,5461 


220. 


.14 


± 


1 


,32 












eg 


01 


:03: 


:19 


.30 


-70: 


:52: 


:26. 


.6 


95, 


93 


164.98 


47874, 


,5479 


214. 


,82 


± 


0, 


.93 




































48607, 


,5636 


215 


.14 


± 





,85 


214.99 


± 


0.63 


0.06 


1 


eh 


01 


:03: 


:21. 


90 


-70: 


:52: 


:04. 


.0 


79, 


,54 


151.74 


47874, 


,5497 


228. 


.24 


± 


0. 


.96 












ej 


01 


:03: 


:23. 


,65 


-70: 


:52: 


:01. 


,4 


81. 


.80 


145.56 


47874, 


,5513 


234, 


.67 


± 


0. 


,77 












ek 


01 


:03: 


:23. 


,77 


-70: 


:51: 


:49. 


.8 


72. 


.92 


140.01 


47874, 


,5530 


224, 


.68 


± 


1, 


37 












em 


01 


:03: 


:24. 


,80 


-70: 


:51: 


:51. 


.5 


77, 


53 


137.95 


47874. 


,5547 


239, 


.52 


± 


0. 


,78 




































48605. 


,6049 


237 


.87 


± 


1. 


.03 


238.92 


± 


0.62 


1.63 


1 



Table 2 (cont.) 
Radial Velocities 
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Star a 5 R PA HJD v r v r \ v Comments 





(2000.0) 


(2000.0) 


(arcsec) 


(°) 


-2400000 


(km 


s 




(km 


s 


!) 






en 


01 


:03: 


25.95 


-70: 


:52:03.8 


90.53 


140.51 


47874.5565 


238.43 


± 


1.11 




























48605.6073 


236.81 


± 


1.12 


237.63 


± 


0.79 


1.06 


1 


ep 


01 


;03: 


28.45 


-70: 


:51:29.4 


78.39 


116.92 


47874.5650 


227.74 


± 


0.96 












eq 


01 


;03: 


25.59 


-70: 


:51:14.3 


59.44 


110.04 


47874.5716 


223.95 


± 


2.08 












er 


01 


;03: 


31.94 


-70: 


:51:15.8 


89.79 


104.12 


47874.5741 


237.84 


± 


2.58 




























48605.6107 


234.44 


± 


1.79 


235.54 


± 


1.47 


1.17 


1 


es 


01 


:03: 


18.36 


-70: 


:51:50.3 


59.88 


160.23 


47874.6101 


217.60 


± 


1.87 












et 


01 


:03: 


17.02 


-70: 


:51:43.9 


51.79 


164.69 


47874.6118 


211.16 


± 


1.53 




























48607.5666 


210.09 


± 


1.34 


210.55 


± 


1.01 


0.28 


1 


eu 


01 


:03: 


13.64 


-70: 


:50:58.7 


5.59 


211.86 


47874.6190 


232.54 


± 


0.98 












ev 


01 


:03: 


12.46 


-70: 


:50:59.8 


10.53 


236.26 


47874.6205 


230.00 


± 


0.88 












ew 


01 


:03: 


10.83 


-70: 


:50:58.5 


17.39 


254.82 


47874.6222 


231.64 


± 


1.15 












f 


01 


:03: 


11.19 


-70: 


:50:46.3 


16.85 


297.00 


47867.5387 


230.67 


± 


0.93 




























47868.5606 


232.80 


± 


1.15 




























47869.5609 


233.23 


± 


0.91 




























47870.5329 


230.69 


± 


1.11 




























47872.5794 


232.42 


± 


1.06 




























47873.6385 


232.65 


± 


1.20 




























47874.6239 


233.12 


± 


0.97 




























48605.5998 


232.95 


± 


1.07 




























48605.6732 


231.12 


± 


1.20 




























48607.5389 


231.67 


± 


1.40 


232.18 


± 


0.34 


8.63 


9 


g 


01 


;03: 


08.88 


-70: 


:50:40.0 


29.84 


297.86 


47867.5405 


214.15 


± 


0.94 












h 


01 


;03: 


10.21 


-70: 


:50:36.3 


26.55 


311.66 


47867.5422 


219.09 


± 


0.94 




























47869.5444 


220.12 


± 


0.75 


219.72 


± 


0.59 


0.73 


1 


j 


01 


:03: 


10.68 


-70: 


:50:32.5 


27.70 


320.75 


47867.5441 


231.75 


± 


0.78 












k 


01 


;03: 


10.46 


-70: 


:50:30.0 


30.33 


322.15 


47867.5459 


225.17 


± 


1.08 












m 


01 


:03: 


11.39 


-70: 


:50:37.6 


21.54 


319.37 


47867.5531 


222.99 


± 


1.07 












n 


01 


:03: 


13.94 


-70: 


:50:40.1 


13.93 


353.91 


47867.5550 


216.43 


± 


1.35 












P 


01 


:03: 


13.60 


-70: 


:50:36.7 


17.54 


349.65 


47867.5568 


228.24 


± 


0.94 




























47873.6292 


229.18 


± 


1.08 


228.65 


± 


0.71 


0.43 


1 


q 


01 


:03: 


13.68 


-70: 


:50:31.8 


22.32 


352.91 


47867.5585 


227.56 


± 


0.80 












r 


01 


:03: 


15.06 


-70: 


:50:32.1 


22.22 


10.4 


47867.5603 


217.66 


± 


0.97 




























47873.6315 


219.93 


± 


1.09 


218.66 


± 


0.72 


2.42 


1 


s 


01 


:03: 


15.99 


-70: 


50:43.5 


13.54 


39.5 


47867.5623 


222.53 


± 


0.80 




























47873.6335 


222.08 


± 


0.87 


222.32 


± 


0.59 


0.14 


1 


t 


01 


:03: 


14.06 


-70: 


:50:51.2 


2.89 


342.15 


47867.5645 


225.73 


± 


1.18 












u 


01 


:03: 


15.19 


-70: 


:51:03.4 


10.54 


153.68 


47867.5664 


232.50 


± 


0.86 












V 


01 


:03: 


15.55 


-70: 


:50:59.6 


8.57 


131.24 


47867.5685 


214.04 


± 


1.31 












w 


01 


:03: 


20.04 


-70: 


:50:54.8 


28.55 


91.7 


47867.5702 


215.61 


± 


0.67 












X 


01 


:03: 


17.41 


-70: 


:50:54.8 


15.62 


93.1 


47868.5255 


227.39 


± 


0.82 




























47873.6353 


226.56 


± 


1.32 


227.16 


± 


0.70 


0.29 


1 


y 


01 


:03: 


17.70 


-70: 


:51:05.4 


20.52 


123.93 


47868.5273 


219.28 


± 


0.71 




























47873.6368 


218.91 


± 


0.81 


219.12 


± 


0.53 


0.12 


1 


z 


01 


:03: 


16.42 


-70: 


:51:07.1 


16.97 


140.80 


47868.5291 


215.14 


± 


1.37 




























47870.5382 


215.51 


± 


0.88 


215.40 


± 


0.74 


0.05 


1 



Table 3 

Expected Fraction of Velocity Differences > (6, 8) km 



x b 




e = 


= 




e = 


/(e) 


o nn 


(0 

v 


DDI 

UUl, 


n ooo^i 


(0 

\ 


DDI 

UUl 


n ooo^i 

U.UUU J 


n o^ 


(0 

V 


099 


o ni 


I'D 

V 


D14 


001 


n i n 
u. ru 








1 U 


09/1 


01 Q"\ 

u.ury j 


0.15 


(o 


071, 


0.063) 


(o 


039 


0.031) 


0.20 


(o 


096, 


0.085) 


(o 


049 


0.040) 


0.25 


(o 


121, 


0.108) 


(o 


069 


0.055) 


0.30 


(o 


147, 


0.130) 


(o 


078 


0.065) 


0.35 


(o 


171, 


0.150) 


(o 


096 


0.079) 


0.40 


(o 


199, 


0.177) 


(o 


107 


0.090) 


0.45 


(o 


221, 


0.195) 


(o 


132 


0.111) 


0.50 


(o 


249, 


0.224) 


(o 


139 


0.118) 







Table 4 








Mass Bins 




oin 




Umax 


— ri 

Contents 




C\Jf \ 

l M oJ 


CM ^ 

l M oJ 




1 


0.15 


0.20 


MS 


2 


0.20 


0.25 


MS 


3 


0.25 


0.30 


MS 


4 


0.30 


0.35 


MS 


5 


0.35 


0.40 


MS 


6 


0.40 


0.45 


MS 


7 


0.45 


0.50 


MS 


8 


0.50 


0.55 


MS & WD 


9 


0.55 


0.60 


MS 


10 


0.60 


0.65 


MS 


11 


0.65 


0.70 


MS & WD 


12 


0.70 


0.75 


MS 


13 


0.75 


0.80 


MS 


14 


0.80 


0.84 


MS 


15 


0.84 


0.87 MS & RG & HB 


16 


0.87 


8.00 


WD 



Notes: MS = main sequence stars 
WD = white dwarf stars 
RG = red giant stars 
HB = horizontal branch stars 



Table 5 

King-Michie - Fitted Parameters - Model A 
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Surface Photometry 



Velocities 



P(IC-ivi) 



x 



Wo 



(pc) 



{y = 22) 



(km s" 1 ) 



ISO 




7.5 ± 


1 


0.56 ±0.02 


10 




7.4 ±0 


1 


0.59 ±0.02 


5 




6.8 ±0 


1 


0.70 ±0.02 


3 




6.0 ±0 


1 


0.80 ±0.03 


ISO 


0.0 


9.1 ± 


2 


0.36 ±0.02 


10 


0.0 


8.7 ±0 


1 


0.39 ±0.02 


5 


0.0 


7.9 ±0 


1 


0.50 ±0.02 


3 


0.0 


6.7 ±0 


1 


0.64 ±0.02 


ISO 


0.5 


10.0 ±0 


2 


0.39 ±0.03 


10 


0.5 


9.6 ±0 


2 


0.43 ±0.02 


5 


0.5 


8.4 ±0 


1 


0.57 ±0.02 


3 


0.5 


7.1 ±0 


1 


0.71 ±0.02 


ISO 


1.0 


11.1 ± 


3 


0.48 ± 0.03 


10 


1.0 


10.9 ±0 


2 


0.52 ±0.02 


5 


1.0 


8.9 ±0 


1 


0.70 ±0.02 


3 


1.0 


7.4 ±0 


1 


0.84 ±0.02 


ISO 


1.5 


12.1 ±0 


3 


0.58 ±0.03 


10 


1.5 


10.9 ±0 


2 


0.70 ±0.02 


5 


1.5 


8.9 ±0 


2 


0.87 ±0.02 


3 


1.5 


7.4 ±0 


1 


0.99 ±0.03 


ISO 


2.0 


13.2 ±0 


3 


0.67 ±0.02 


10 


2.0 


10.7 ±0 


3 


0.88 ±0.03 






60 


0.91 





51 


0.96 


1 


02 


0.37 


2 


28 


0.00 


1 


83 


0.01 





77 


0.74 





60 


0.92 


1 


38 


0.08 



205.69 


0.02 


205.99 


0.07 


206 55 


26 


207.27 


0.72 


206.09 


0.09 


206.69 


0.22 


207.02 


0.59 


207.47 


0.88 


205.54 


0.01 


206.20 


0.13 


206.63 


0.34 


207.09 


0.66 



2.0 
2.0 



8.8 ±0.2 
7.2 ±0.1 



1.02 ±0.03 
1.10 ±0.03 



49 ± 4 
76 ±13 

148 

133 



+67 
-51 
+85 
-47 



56 ± 5 
87 ±14 
217 + 70 

ZLI -20 



173 



+63 
-47 



57 ± 5 
99 ±20 
192 tlo 
159 ±g 

59 ± 5 
212 ±60 
156 
140 



+40 
-60 
+61 
-46 



72 ± 6 
151 



+30 
-50 

145 + 60 



128 



-30 
+56 
-26 



105 ± 10 

136 ± 40 

137 + 40 

io ' -20 

68 ±20 



1.32 
0.64 
0.93 
2.06 



0.12 
0.86 
0.49 
0.00 



8.9 ±0.5 
9.2 ±0.5 
9.9 ±0.5 
11.2 ±0.6 

8.5 ±0.4 
9.2 ±0.5 
9.9 ±0.5 
11.1 ±0.6 

7.7 ±0.4 
8.4 ±0.4 
9.2 ±0.5 
10.4 ±0.6 



0.66 


0.85 


7.0 ±0.4 


204.99 


0.00 


0.52 


0.96 


7.5 ±0.4 


205.39 


0.01 


2.02 


0.00 


8.7 ±0.4 


206.08 


0.10 


3.51 


0.00 


9.9 ±0.5 


206.59 


0.36 



0.60 


0.91 


1.28 


0.16 


3.79 


0.00 


5.55 


0.00 


0.99 


0.41 


2.98 


0.00 


5.64 


0.00 


7.68 


0.00 



6.7 ±0.3 

7.4 ±0.4 

8.5 ±0.4 

9.8 ±0.5 

6.7 ±0.3 

7.5 ±0.4 

8.6 ±0.4 
10.1 ±0.5 



204.74 


0.00 


205.01 


0.00 


205.55 


0.02 


206.12 


0.15 


204.70 


0.00 


204.83 


0.00 


205.26 


0.00 


206.08 


0.11 



Table 6 

King-Michie - Derived Parameters - Model A 
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(*») 


X 


(w«pc- 3 ) 


Ly 
10 5 L VQ 


Po 

(m pc" 3 ) 


M 

(1O 5 M ) 


Population 
(M/L v )o (M/L v ) 
(M/L v )o 


Dynamic 
(M/L v ) (M/L v ) 
(M/Mq 


ISO 




3.0 ± 0.2 


1.66 ± 0.05 


4.2 ± 0.5 


2.3 ± 0.3 






1.4 ± 0.2 


1.4 ± 0.2 


10 




2.8 ±0.2 


1.72 ±0.06 


4.1 ±0.5 


2.5 ±0.3 






1.4 ± 0.2 


1.4 ± 0.2 


5 




2.2 ±0.2 


1.73 ±0.06 


3.3 ±0.4 


2.7 ±0.3 






1.6 ±0.2 


1.6 ±0.2 


3 




1.9 ± 0.2 


1.65 ± 0.05 


3.3 ± 0.5 


2.9 ± 0.3 






1.7 ±0.2 


1.7 ±0.2 


ISO 


0.0 


3.7 ± 0.4 


1.57 ± 0.04 


9.4 ± 1.5 


2.3 ± 0.2 


7.03 


4.21 


or i r\ o 

2.5 ± 0.3 


1.5 ± 0.2 


10 


0.0 


3.5 ±0.3 


1.65 ±0.05 


9.0 ± 1.0 


2.6 ±0.3 


6.88 


4.21 


2.6 ± 0.3 


1.6 ± 0.2 


5 


0.0 


2.6 ±0.3 


1.77 ±0.05 


6.4 ± 1.0 


3.0 ±0.3 


6.03 


4.21 


2.5 ±0.3 


1.7 ±0.2 


3 


0.0 


2.1 ± 0.2 


1.71 ± 0.05 


5.1 ± 0.6 


3.2 ± 0.4 


5.44 


4.21 


2.4 ±0.3 


1.9 ±0.2 


ISO 


0.5 


3.6 ± 0.3 


1.60 ± 0.04 


6.3 ± 1.0 


2.6 ± 0.3 


3.45 


3.17 


1.8 ± 0.2 


1.6 ± 0.2 


10 


0.5 


3.2 ±0.2 


1.67 ±0.05 


6.2 ±1.0 


3.0 ±0.3 


3.38 


3.17 


1.9 ± 0.2 


1.8 ± 0.2 


5 


0.5 


2.3 ±0.2 


1.72 ±0.05 


4.4 ±0.5 


3.4 ±0.4 


2.98 


3.17 


1.9 ±0.2 


2.0 ±0.2 


3 


0.5 


1.9 ± 0.2 


1.67 ± 0.05 


3.6 ± 0.4 


3.5 ± 0.4 


2.77 


3.17 


1.9 ±0.2 


2.1 ±0.3 


ISO 


1.0 


3.1 ± 0.3 


1.65 ± 0.05 


3.6 ± 0.5 


3.2 ± 0.4 


1.78 


3.07 


1.1 ± 0.1 


2.0 ± 0.2 


10 


1.0 


2.6 ±0.2 


1.77 ±0.06 


3.4 ±0.4 


4.2 ±0.5 


1.66 


3.07 


1.3 ± 0.1 


2.4 ± 0.3 


5 


1.0 


2.0 ±0.1 


1.65 ±0.05 


2.6 ±0.3 


4.0 ±0.5 


1.63 


3.07 


1.3 ± 0.1 


2.4 ± 0.3 


3 


1.0 


1.7 ± 0.1 


1.58 ± 0.04 


2.4 ± 0.3 


4.0 ± 0.5 


1.62 


3.07 


1.4 ±0.2 


2.6 ±0.3 


ISO 


1.5 


2.8 ± 0.2 


1.71 ± 0.05 


2.2 ± 0.3 


4.6 ± 0.5 


1.06 


3.50 


0.8 ±0.1 


2.7 ±0.3 


10 


1.5 


2.1 ± 0.1 


1 70 ± 05 


1.8 ± 0.2 


5.1 ± 0.7 


1.01 


3.50 


0.9 ±0.1 


3.0 ±0.4 


5 


1.5 


1.7 ±0.1 


1.56 ±0.04 


1.6 ±0.2 


4.8 ±0.6 


1.07 


3.50 


1.0 ±0.1 


3.1 ±0.4 


3 


1.5 


1.5 ±0.1 


1.49 ± 0.04 


1.6 ±0.2 


4.6 ±0.6 


1.16 


3.50 


1.1 ±0.1 


3.1 ±0.4 


ISO 


2.0 


2.4 ±0.2 


1.83 ±0.06 


1.6 ±0.2 


7.8 ± 1.0 


0.70 


4.50 


0.7 ±0.1 


4.3 ±0.5 


10 


2.0 


1.7 ±0.1 


1.62 ±0.05 


1.2 ±0.1 


6.3 ±0.5 


0.79 


4.50 


0.7 ±0.1 


3.9 ±0.5 


5 


2.0 


1.5 ±0.1 


1.47 ±0.04 


1.2 ±0.1 


5.6 ±0.6 


0.92 


4.50 


0.8 ±0.1 


3.8 ±0.4 


3 


2.0 


1.6 ±0.1 


1.45 ±0.04 


1.4 ±0.2 


4.9 ±0.6 


1.14 


4.50 


0.9 ±0.1 


3.4 ±0.4 
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Surface Photometry 



Velocities 



P(IC-ivi) 



x 



Wo 



(pc) 



{y = 22) 



(km s" 1 ) 



206.01 


0.08 


206.61 


0.28 


206.96 


0.55 


207.44 


0.86 


205.48 


0.01 


206.10 


0.10 


206.58 


0.31 


207.04 


0.62 


204.97 


0.00 


205.37 


0.00 


206.03 


0.09 


206.55 


0.34 


204.73 


0.00 


204.99 


0.00 


205.52 


0.02 


206.12 


0.15 


204.70 


0.00 


204.83 


0.00 


205.25 


0.00 


205.82 


0.06 



0.0 
0.0 
0.0 
0.0 



ISO 


1.0 


10 


1.0 


5 


1.0 


3 


1.0 


ISO 


1.5 


10 


1.5 


5 


1.5 


3 


1.5 


ISO 


2.0 


10 


2.0 


5 


2.0 


3 


2.0 



9.1 ±0.2 
8.7 ±0.1 
7.9 ±0.1 
6.7 ±0.1 



0.37 ±0.02 
0.40 ±0.02 
0.52 ±0.02 
0.65 ±0.02 



57 ± 5 

88 ±14 

221 + 70 
zzi _ 20 

169 tf 7 



10.9 ±0.3 
10.8 ±0.2 

8.8 ±0.1 
7.4 ±0.1 

12.1 ±0.3 
10.8 ±0.2 

8.9 ±0.2 
7.4 ±0.1 

13.3 ±0.3 
10.6 ±0.3 
8.8 ±0.2 
7.4 ±0.1 



0.50 ±0.03 
0.54 ±0.02 
0.72 ± 0.02 
0.85 ±0.02 

0.58 ±0.03 
0.72 ±0.02 
0.89 ±0.02 
1.00 ±0.03 

0.67 ±0.02 
0.88 ±0.03 
1.03 ±0.03 
1.14 ±0.03 



1.70 
0.73 
0.63 
1.48 



0.02 
0.79 
0.90 
0.05 



ISO 





5 


10.0 ±0 


2 


0.40 ±0.03 


58 ± 5 


1 


20 


0.19 


10 





5 


9.6 ±0 


2 


0.45 ±0.02 


101 ±20 





61 


0.90 


5 





5 


8.4 ±0 


1 


0.58 ±0.02 


184 t£ 


1 


03 


0.36 


3 





5 


7.1 ±0 


1 


0.72 ±0.02 


157 tit 


2 


22 


0.00 



59 ± 5 





58 


0.93 


195 ± 60 





56 


0.94 


156 tf 


2 


21 


0.00 


145 + _% 


3 


73 


0.00 


75 ± 6 





60 


0.91 


156 tf 


1 


42 


0.08 


151 tf 


4 


10 


0.00 


-| r)7 +56 
LZI -26 


5 


93 


0.00 


111 ± 10 





94 


0.48 


129 ±40 


3 


35 


0.00 


138 tf 


6 


16 


0.00 


135 tf 


8 


30 


0.00 



8.4 ±0.4 

9.1 ±0.5 

9.8 ±0.5 
11.0 ±0.6 

7.6 ±0.4 
8.3 ±0.4 

9.2 ±0.5 
10.4 ±0.6 

7.0 ±0.4 

7.5 ±0.4 

8.6 ±0.4 

9.9 ±0.5 

6.7 ±0.3 

7.3 ±0.4 

8.4 ±0.4 

9.7 ±0.5 

6.6 ±0.3 

7.5 ±0.4 

8.6 ±0.4 

9.8 ±0.5 
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(*») 


X 


(w«pc- 3 ) 


Ly 
10 5 L VQ 


Po 

(m q pc" 3 ) 


M 

(1O 5 M ) 


Population 
(M/L v )o (M/L v ) 
(M/L v )o 


Dynamic 
(M/L v ) (M/L v ) 
(M/L v )o 


ISO 


0.0 


3.7 ± 0.4 


1.58 ± 0.04 


8.8 ± 1.5 


2.4 ± 0.2 


4.06 


2.56 


2.4 ± 0.3 


1.5 ± 0.2 


10 


0.0 


3.5 ±0.3 


1.66 ±0.05 


8.5 ±1.0 


2.6 ±0.3 


3.97 


2.56 


2.4 ± 0.3 


1.6 ± 0.2 


5 


0.0 


2.6 ±0.3 


1.77 ±0.05 


6.0 ±1.0 


3.0 ±0.3 


3.47 


2.56 


2.4 ±0.3 


1.7 ±0.2 


3 


0.0 


2.1 ± 0.2 


1.70 ± 0.05 


4.8 ± 0.6 


3.2 ± 0.4 


3.16 


2.56 


2.3 ±0.3 


1.9 ±0.2 


ISO 


0.5 


3.6 ± 0.3 


1.61 ± 0.04 


6.0 ± 1.0 


2.6 ± 0.3 


2.01 


1.96 


1.7 ± 0.2 


1.6 ± 0.2 


10 


0.5 


3.2 ±0.2 


1.68 ±0.05 


5.7 ± 1.0 


3.0 ±0.3 


1.96 


1.96 


1.8 ± 0.2 


1.8 ± 0.2 


5 


0.5 


2.3 ±0.2 


1.71 ±0.05 


4.1 ±0.5 


3.4 ±0.4 


1.75 


1.96 


1.8 ±0.2 


2.0 ±0.2 


3 


0.5 


1.9 ± 0.2 


1.66 ± 0.05 


3.4 ± 0.4 


3.5 ± 0.4 


1.64 


1.96 


1.8 ±0.2 


2.1 ±0.3 


ISO 


1.0 


3.0 ± 0.3 


1.66 ± 0.05 


3.3 ± 0.5 


3.3 ± 0.4 


1.04 


1.91 


1.1 ± 0.1 


2.0 ± 0.2 


10 


1.0 


2.6 ±0.2 


1.75 ±0.06 


3.2 ±0.4 


4.2 ±0.5 


1.00 


1.91 


1.2 ± 0.1 


2.4 ± 0.3 


5 


1.0 


1.9 ±0.1 


1.64 ±0.05 


2.4 ±0.3 


4.0 ±0.5 


0.97 


1.91 


1.3 ± 0.1 


2.4 ± 0.3 


3 


1.0 


1.7 ± 0.1 


1.57 ± 0.04 


2.2 ± 0.3 


4.0 ± 0.5 


0.97 


1.91 


1.3 ±0.2 


2.6 ±0.3 


ISO 


1.5 


2.8 ± 0.2 


1.72 ± 0.05 


2.2 ± 0.3 


4.7 ± 0.5 


0.65 


2.22 


0.8 ±0.1 


2.7 ±0.3 


10 


1 5 


2.0 ± 0.1 


1 69 ± 05 


1 7 ± 2 


5 1 ± 7 


0.61 


2.22 


0.9 ±0.1 


3.0 ±0.4 


5 


1.5 


1.6 ±0.1 


1.54 ±0.04 


1.5 ±0.2 


4.8 ±0.6 


0.65 


2.22 


0.9 ±0.1 


3.1 ±0.4 


3 


1.5 


1.5 ±0.1 


1.47 ±0.04 


1.6 ±0.2 


4.6 ±0.6 


0.73 


2.22 


1.0 ±0.1 


3.1 ±0.4 


ISO 


2.0 


2.4 ±0.2 


1.82 ±0.06 


1.6 ±0.2 


8.1 ± 1.0 


0.43 


2.90 


0.7 ±0.1 


4.5 ±0.5 


10 


2.0 


1.7 ±0.1 


1.59 ±0.05 


1.2 ±0.1 


6.2 ±0.5 


0.50 


2.90 


0.7 ±0.1 


3.9 ±0.5 


5 


2.0 


1.5 ±0.1 


1.44 ± 0.04 


1.2 ±0.1 


5.5 ±0.6 


0.57 


2.90 


0.8 ±0.1 


3.8 ±0.4 


3 


2.0 


1.4 ±0.1 


1.36 ±0.04 


1.2 ±0.2 


5.3 ±0.6 


0.67 


2.90 


0.9 ±0.1 


3.9 ±0.4 



